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Abstract 



A generalization of the Bethe ansatz equations is studied, where a scalar 
two-particle S-matrix has several zeroes and poles in the complex plane, 
as opposed to the ordinary single pole/zero case. For the repulsive case 
(no complex roots), the main result is the enumeration of all distinct 
solutions to the Bethe equations in terms of the Fuss-Catalan numbers. 
Two new combinatorial interpretations of the Fuss-Catalan and related 
numbers are obtained. On the one hand, they count regular orbits of the 
permutation group in certain factor modules over Z M , and on the other 
hand, they count integer points in certain M-dimensional polytopes. 
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Introduction 



The modern techniques of quantum integrability rely on what is known today as the Bethe 
ansatz, a powerful tool, first introduced by H. Bethe |T] and refined by many authors 
throughout the years, see the summarising monographs [2j [3j Hj. The idea behind the 
Bethe ansatz is to reduce, in one way or another, the spectral problem for the Hamiltonian 
in the M-particle sector to the one of solving a set of M simultaneous equations (algebraic 
or transcendental, depending on the model) for particles' rapidities v n . 

In the case when the particles carry no internal degrees of freedom (are spinless), those 
equations, referred to as Bethe equations, can be cast in the form 

M 

F ( v n) IJ S(v n -v m ) = 1, n = l,...,M, (0.1) 

m=l 

and interpreted as the periodicity conditions for the wave function of the system of M 
one-dimensional particles put in a periodic box, in the assumption that all multi-particle 
collisions can be reduced to a sequence of two-particle collisions (factorisation of S -matrix 
hypothesis). The function F can be interpreted then as a one-particle travel phase factor, 
and S two-particle scattering phase factor. 



Most of the solutions of the equations (0.1), referred to in the present paper as weak 



generalized Bethe equations (weak GBE), correspond to the discrete eigenstates of the 
integrable system. To eliminate the nonphysical solutions one needs, as explained in 
section [TJ to supplement (0.1) with additional conditions. Such corrected GBE will be 



called strong GBE. 

Although F can take a wide range of forms depending on the chosen integrable model, 
the expression for S is much less variable, being restricted by the quantum group governing 
the integrability of the model. In the most of the presently known examples, the S-factor 
has only a single pole and a single zero in the complex plane (modulo periodicity for 
trigonometric or elliptic functions). The only exception known so far is the Izergin- 
Korepin model [HIE]. 

In the present paper we start a systematic study of the Bethe equations, where the 
function S has K poles and K zeroes. Our primary interest is the case K > 1, referred to 
as the generalized Bethe equations (GBE), as opposed to the well-studied ordinary case 
K — 1. Our main motivation comes from the recent progress on AdS/CFT correspondence 
in string theory that produced a similar kind of Bethe equations |S]- It has to be 
emphasized that in contrast to the "asymptotic Bethe ansatz" approach as in [Zl IS], 



where the equations (0.1 ) hold in the infinite volume limit only, we assume here that they 
are satisfied exactly in the finite volume as well. 

We shall concentrate on the simplest task: counting the number of solutions to the 
strong Bethe equations. Even in the ordinary case of a single pole/zero 5-function, the 
combinatorics of Bethe ansatz related to the so-called "string hypothesis" and "complete- 
ness problem" has been found to be rewardingly rich and interesting [91 [10] . Mathemat- 
ically, the task is to count solutions of a special system of algebraic equations. We hope 
that our results present interest not only for the specialists in integrable systems but also 
in combinatorics. 
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To avoid complications caused by complex solutions we focus on the repulsive case, 
when all roots v n are real. In this case, the strong conditions are equivalent simply to 
discarding multiple roots of the Bethe equations. We consider two cases: the functions F 
and S being ratios of trigonometric or ordinary polynomials of degree N for F and degree 
K for S. The number of distinct solutions to the strong GBE is then denoted C^' N \ 

To count the solutions of GBE we use the fact that they correspond to the minima of 
certain convex function (Yang- Yang function). In the trigonometric case, we establish a 
one-to-one correspondence between the set of distinct solutions to the GBE and the set 
of regular orbits of the permutation group &m in certain factor modules over Z M . In the 
rational case, we bijectively label the solutions with the integer sequences I\ < ■ ■ ■ < Im 
satisfying certain linear inequalities. For a more detailed description see respectively 
theorems 1 and 2 below. 

For the numbers C^' N ^ that enumerate both combinatorial objects and therefore the 
solutions of strong GBE as well, we then obtain an explicit formula 

i(k,n) N [KM + N-1\ N [KM + K 



° M ~M\ M-l )~KM + n{ M ) (0 ' 2) 

that constitues the main result of our paper. 

The respective generating function turns out to be the N th power of the generating 
function for the Fuss-Catalan numbers [TTJ [T2] corresponding to C^f i n our notation: 

oo , N 

C^ N \t) = ]T C[f N) t M = lc^\t) J . (0.3) 

M=0 ^ ' 

The particular case C M ^' corresponds to the famous Catalan numbers [13] that are 
known to have over 120 different combinatorial interpretations [T^j . 

Two theorems formulated below provide two new combinatorial interpretations of the 
numbers C^' , including Fuss-Catalan numbers C^' 1 ^ and Catalan numbers C^' 1 ^ as 
particular cases, obtained as a result of our study of GBE. 

Theorem 1. Let N, M and K be some positive integers. Two M -tuples of integers 
(J 1; . . .Im) ore considered equivalent if one can be mapped into another by a finite com- 
position of basic transformations 7„ (n = 1, . . . , M) 

ln :I n v->I n + N+(M- l)K and 7n : I m \-> I m - K m^n, (0.4) 

An equivalence class is called regular if it contains only M -tuples of pairwise distinct 
integers (Ii, . . . Im)- 

Then the number of all regular equivalence classes is given by Ml C^ ,N \ 

Theorem 2. Let N , M and K be some positive integers, then the number of different 
choices of strictly increasing integer sequences I\ < ■ ■ ■ < Im satisfying to 2M constraints 

rij ~ - K < h + I 2 + ... + I r , (0.5a) 

I M -r+i + ... + Im-i + Im < rMK- ^ + ^ K + rN-l, (0.5b) 
for r — 1, . . . , M equals C^' N ^ . 
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The paper is organized as follows. In Section [TJ we give a precise description of the 
class of functions F and S that we consider. We then define the strong generalized Bethe 
equations (GBE) and prove certain elementary properties of the solutions to the weak and 
strong GBE associated with our choice of F and S. In Section [2j we reformulate GBE 
as the minimality condition for a convex Yang- Yang function and establish a one-to-one 
correspondence between the set of distinct solutions to the strong GBE and the sets of 
integers described in theorems [T] and [2] Then, we distinguish between the trigonometric 
and rational cases. In Section |3j we apply the theory of commutative rings to count 
the distinct solutions to the strong GBE in the trigonometric case and prove theorem 
[TJ Finally, in Section [4], we analyze the limit when the period of trigonometric functions 
goes to infinity and show that the number of solutions in the rational case is the same as 
in the trigonometric case. Theorem [2] is then obtained as a corollary of the established 
bijection. 

The important question as to what kind of integrable models and what quantum 
groups might correspond to our generalized Bethe equations is left unresolved. We plan 
to address it in our subsequent publications. 



1 The general setting 

Throughout the paper, we consider two cases of interest. We assume that F and S are 
either realized as ratios of ordinary polynomials or as ratios of trigonometric polynomials. 
We will refer to the first case as the rational case and to the second one as the trigonometric 
case. 

Let ip(u) = u in the rational case and ip{u) = sin(w) in the trigonometric one. The 
functions F and S are characterized by natural numbers N, K and complex parameters 
Foc,£n, Vj C and are defined as 

N 

F(u) = F x f[ y' U + l n) 1^1 = 1, lm(£„)>0, ne{l,...,N}, (1.1a) 
= II TTX^Y ' Mvj)>0, je{l,...,K}. (l.lb) 

fJl ^( u + Vj) 

Here z stands for the complex conjugate of z. The parameters rjj are chosen in such 
a way that the sets of zeroes and poles of S are both symmetric with respect to the 
imaginary axis: 

{v} = {~v}- (1.2) 
To avoid degeneracy of solutions to GBE, we demand that F^ ^ 1, or, equivalently, 
Foo = e -2 ' 7 ^ 00 , with — 1 < (foe < 0. Later, in some cases we impose more restrictive 
conditions on 

Finally, in the trigonometric case, we also assume that the parameters £ n and rjj are 
located in the strip —ir/2 < Re(z) < 7r/2. 

The above restrictions on F and S are chosen in such a way that the phase factors, 
both trigonometric and rational, have the following natural, for a physicist, properties 



that are easy to derive from (1.1) and (1.2). 
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1. Unitarity. 



F(u) = F{u)-\ S(u) = S{u)-\ u e C. (1.3) 

2. Parity. S respects the space-, or P-parity: 

S(-u) = S{u)-\ (1.4a) 

and the time-, or T-parity: 

S{uj = S(-u). (1.4b) 

3. Repulsion condition: 

\F(u)\ > 1, \S{u)\>l, lmw>0, (1.5a) 
\F(u)\ < 1, \S{u)\<l, lmw<0. (1.5b) 

Remark. The unitarity condition implies that |P(w)| = 1 = ^(w)! for uel Also 
it is easy to see that the two parity conditions (1.4a) and (1.4b) are equivalent as long as 
one assumes unitarity. 

It is a well known fact in the literature on Bethe ansatz that the Bethe equations (0.1 ), 
taken literally as they are written, have too many "false" solutions that do not correspond 
to any physical states. A common trick to eliminate them from consideration is to replace 
(0.1) with a different, more restrictive problem for v m in terms of the Baxter polynomial: 

M 

Q(u\ v) = Y[ ip(u - v m ), v = (vi, . . . ,v M ). (1.6) 

m=l 

The problem for Q(u \ v) described below generalizes Baxter's TQ equation known in 
the ordinary K = 1 case [2]. An important difference is that for K > 1 the equation for 
Q is not linear. 

Let z = (zi, . . . , Zm) £ C A/ and Q{r ? },{5}('^| z) be defined as 

K K 

Q {r , m} (u\z) = a^\{0)l[Q(u + r ]k \z) + S(u\ {£}) J] Q(u + r] k \ z) . (1.7) 

k=l k=l 

with 

N N 

a(u\ {£}) = Poo J] ^ u + £n) , 8(u\ {£}) = ("1) A ' +1 II ^ u + O • ( L8 ) 

n—1 n=l 



We define the strong GBE to be the set of M conditions on the coordinates of the 
vector v = (v\, . . . ,v m ) necessary and sufficient that the function 

M ^ aw.«>W (L9) 

Q(u\ v) 

is entire. 

One can see that the above definition implies that Q{??},{£}(w| v) vanishes at u — Vj 
and hence the roots Vj always satisfy the weak GBE (0.1). However, the definition of the 
strong GBE is slightly more restrictive as it provides additional equations in the case of 
multiple roots [ie Vj = v k for some j ^ k). 
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1.1 General properties of solutions 

In this subsection, we discuss general properties of the solutions to the weak and strong 



GBE with F and S given by (1.1). We first show that any solution to the weak, and 
hence to the strong GBE as well, has necessarily real roots Vj. Then, we show that 
for the strong GBE there are no "double" roots (ie Vj ^ Vk for j ^ k). For the both 
properties the repulsion condition (1.5) is crucial. The proofs follow closely the standard 
ones [3j H] for the case K = 1. Note that, until section [3} we ignore the 7r-periodicity in 



the trigonometric case and assume that Vj G K. If in a pair of solutions v and v' we have 
v'„ = v n + Ti the solutions are treated as different ones. 



Proposition 1. Let F and S satisfy the repulsion conditions (|l.5|). Then all of the 
coordinates v n of a solution v 



vm) to the weak GBE (0.1) are real. 



Proof. The proof goes as in [3J. One orders the roots v n with respect to their increasing 
imaginary part: 

Imt>i < Imt>2 < • • • < Im%. (1-10) 
Then lm(vM — v m ) > Vm and therefore \S(vm — v m )\ > 1 due to the repulsion 



condition (1.5). The GBE (0.1) then imply \F(vm)\ < 1 that is Im% < 0, as F satisfies 
the repulsion condition. Similarly, setting n = 1, we obtain that Im(vi — v m ) < Vm 
leading to \S(vi — v m )\ < 1 and therefore |F(ui)| > 1, that is to say Imt>! > 0. 

The two inequalities Imi^ < and Imfi > along with (1.10) lead to lmv n = for 
all n. m 



Proposition 2. [The exclusion principle] Let F, S be as in (1.1), and v be a solution to 
the strong GBE. Then the coordinates of v are pairwise distinct, ie Vj ^ Vk for j ^ k. 

Proof. Again, we follow closely a proof from [3J for the case K = 1. Assume that there 
exists a solution v = (v-i, . . . , Vm) to the strong GBE with a double root Vk = Vj for some 



j ^ k. The condition for the expression (1.9) to be entire implies that u = v k is a double 



root of Q{r]},{s,}(u\ v) (1.7) and thus 

d 

du 



0. 



(1.11) 



u=v k 



By proposition [TJ all the roots Vk are real. Therefore, 



K M 



a 



i=\ m=l 



and one can apply the weak GBE (0.1) to recast the derivative of Q{ri},{£} as 



du 



K M 



u=v k 



=1 m=l 



X 



N 

E 



K M 



ip(2i Im^) 



=1 m- 



^ \ifj(v k -v m + r)e)\' 



. (1.12) 
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The factor appearing in the second line is purely imaginary with a strictly positive 
imaginary part and the pre-factor does not vanish. This contradicts the assumption that 

&{vHt}( v * I v ) = °- ■ 
We have thus proven that in the repulsive case the strong GBE are equivalent to the 



weak GBE (0.1) supplemented with the condition Vj ^ V/- for j ^ k. From now on, we 



can focus on the analysis of the weak GBE (0.1) only, simply discarding the solutions 
with multiple roots. 

From the physical point of view, the solutions (vi, . . . , vm) differing by a permutation 
of roots Vj are indistinguishable, defining the same physical state. In the repulsive case, 
when all roots are real and distinct, it is sufficent then to consider only the ordered 
solutions that satisfy 

vi<v 2 < ... < v M . (1.13) 



In this text, however, we do not assume condition (1.13), unless stated otherwise. The 



solutions satisfying (1.13) are always explicitely referred to as ordered solutions. 



1.2 Logarithmic Bethe equations 

The logarithmic Bethe equations P, H] provide a way of rewriting the Bethe equations in 
such a way that a given solution v = (vi, . . . , %) becomes labelled by a set of M integers 
I\, . . . , Im- These equations (logBE) are obtained by taking the logarithm of (0.1 ): 



M 



+ ^2 d ( Vn ~ Vm ^ = 27rL " 



n 



1 . . . M , 



m=l 



with ip and 9 related to the logarithms of F and S 

(p(u) = i In F(u) and 9(u) = ilnS'(u) 



;i.i4) 



;i.i5) 



We fix the branch of the logarithm of S(u) by demanding that 8(0) = irK, which is 



possible due to the symmetry ( 1.2 ). As for the branch of In F(u), we just assume that it is 



fixed in some arbitrary way. We also stress that <p and 9 are strictly increasing functions 



<P W 



—i 



N 

E 



ip(2i Img w ) 



K 



>0 and 9'{u) = -iY^ 



> 0. 



;i.i6) 



Moreover, due to (1.2), 9' is even. 



To be able to use the above correspondence of the solutions to GBE and M-tuples 
of integers we need to describe the class of admissible integers I n that can arise in the 



right-hand side of (1.14) and to establish a bijection between those M-tuples and the 
solutions to GBE. 

The construction is different in the trigonometric and rational cases and is described 
in the next section. In the rest of the current section, we present the properties of the 
integers I n that are shared by the trigonometric and rational cases. 
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Proposition3. [Monotonicity principle] . Under the repulsion conditions (1.5), a growing 
sequence of roots v n of the logBE (1.14) corresponds to a growing sequence of integers I n , 
and vice versa: 



V k > Vj 



I, > L. 



;i.i7) 



Proof. Suppose that v 2 > V\. From the strict monotonicty (1.16) of (p(u), resp. 9(u), it 
follows that (p(v2) > <p(vi), resp. 9(v 2 — V\) > 6{v\ — v 2 ) and 9{v 2 — v m ) > 9(v\ — v m ) for 
m > 2. Adding all the inequalities together we obtain I 2 > I±. 

Now assume that I 2 > I\. If V\ = v 2 then clearly I\ = I 2 , that contradicts the 
assumption. Hence v\ ^ v 2 . Also, the variant v 2 < v\ is impossible as, due to the above 
paragraph, it would imply that I\> 1%. Therefore v<i > V\. m 

Proposition 4. Let I = (Ii, . . . ,1m) be M integers such that there exists a solution 
v = (%,..., Vm) to the logBE (1.14). Then v is the only solution corresponding to the 
same I . 



Proof. Note that the logBE (1.14) appear as the extremum condition w n = dW/dv n = 
for the potential W(v) on IR A ^ 

M . M AI M 

W(v) = '$2$(v n ) + - Q(v n -v m )-27rJ2l n v n + 7r(M-l)Kj2v n , (1.18) 



n=l 



m,n=l 



n=l 



n=l 



where $ and are the integrals of ip and 9 (1.15) 

f'U 

$(n) = / ip(t) dt, Q(u) 



9(t) dt . 



;i.i9) 



The calculation of w n uses the identity 9(u) + 9(—u) = 29(0) = 27rK that follows from 
9' being an even function. 

To prove the uniqueness of solutions to ( 1.14 ) for a given choice of J, it suffices to show 
that the potential is strictly convex, and therefore, if it admits a minimum, the latter is 
unique. To prove the strict convexity of W it is enough to check that the Hessian of W 



d 2 W 
dv n dv p 



M 



- e'(v r . 



1.20) 



m=l 



is a positively definite matrix. Evaluating the quadratic form of the Hessian (1.20) on an 
arbitrary nonzero vector g G 1R A/ we obtain 



M 

E 

n,p=l 



d 2 W 
dv n dv. 



M 



M 



M 



9n9 P = + S 6 '( Vn ~ Vm ^ / 6 ' ( Vn ~ V P^ 9n9 P 



n=l 
M 



m=l 
M 



n,p=l 



M 



^2<P f (v n )g* + Y e '( v n-v m )(g n - g m ) 2 > ^2<p'(v n )gfc > 0. (1.21) 



n=l 



71=1 



0. 



Here we used that (p' > and 9' > on K, and that 6 1 ' is even. ■ 
In the K = 1 case the potential W is of course the well-known Yang- Yang function 



S 



2 Solvability of the generalized Bethe equations 



In this section we derive the necessary and sufficient conditions on the integers Ij for 
the corresponding GBE to have a (unique, by proposition [1]) solution. The problem of 
counting of the solutions to GBE will be thus reduced to the enumeration of the allowed 
sets of integers. The description of those sets is different in the trigonometric and rational 
cases. 

2.1 Trigonometric case 

It is useful to introduce an auxiliary function: 

u 

X (u,c) = iln Shl (; + iC l = [ . , Sinh(2c ; -dv, (2.1) 



smitt — ic) J smiu + ic) sin(i> — ic) 



2 



fixing the branch of logarithm by the condition x(— 7r/2) = 0. Note that u \- > x(w, c) is 7r 
quasi-periodic: 

x(u + iik,c) = 2irk + x(u,c), fceZ and such that x(~ u i c ) + x( u i c ) — 2tt . (2.2) 

Moreover, x( u -> c ) grows monotonously on the real axis and its integral has the asymptotic 
behavior 

I x(t,c)dt = u 2 + 0(u), u^±oo. (2.3) 
Jo 

The functions x{ u i c ) allow us to provide a more explicit representation for tp and 9: 



N 



ip{u) = 27^ + X (u + Re(Cn), Im(a)) , where = e" 2 ^, (2.4a) 



n=l 

K 



0(u) = Y^x{u + Re(^), Im(^)). (2.4b) 

3=1 



One can check using (2.2) and (1.2) that indeed one has ^(0) = 7r.fr. 



We are now in position to prove the existence of solutions in the trigonometric case. 

Proposition 5. Let I = (Ji, . . . , Im) € ^ M be any sequence of M integers. Then, in the 
trigonometric case where ip{u) = sin(-u), there exists a unique solution to the logBE. 

Proof. We have already seen in proposition [4] that, if a solution exists, then it is unique as 



it provides the minimum for a strictly convex potential W (1.18). To prove the existence 
of solutions, it is thus enough to show that W grows at infinity. Indeed, W is then 
bounded from below and hence has a local minimum. 



It follows from (2.3) that $ and G, the integrals of ip and 6, have the asymptotic 
behavior 



$( n ) = / v {t) dt = Nu 2 +0{u) , G(u) = / 6{t) dt = Ku 2 +0(u) for u -> ±oo . 

(2.5) 



o 
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Hence, it follows that W(v) has the asymptotics 

M r , M 

K 



W{v) = nY^vI + — K-^ m ) 2 + 0(max|^|). 



(2.6) 



n=l 



m,n=l 



The potential W thus grows in every direction at infinity, and hence W has a local 



minimum. The existence of solutions to (1.14) is thus proven. 



So far we ignored completely the periodicity of trigonometric functions. However, 
the periodicity F(u + ir) = F(u), S(u + n) = S(u) means that two different solutions 
v,v' G IR of the logarithmic Bethe equations (1.14) are equivalent as solutions of the 



weak GBE (0.1) if v — v' G (7rZ) M . To obtain a characterization of the set of distinct 
solutions one needs to express the equivalence of the vectors v modulo the lattice (7rZ) M 
in terms of the integers I. The answer is given below. 
Denote 

a = N+ (M — 1)K, b=-K. (2.7) 



Proposition 6. Let two solutions v, v' G 
single root v n by it: 

v' m =V m + Tfd'rnn, 171=1 

Then for the corresponding integer vectors we have 



M of the logBE fll~l~4| differ by the shift of a 
,M. (2.8) 



I' m = I m + a5 mn + b(l - 5 n 



m 



1 , M. 



(2.9) 



Proof. The statement follows from the logBE (1.14) and from the quasi-periodicity 
conditions 

ip(u + n) = 2nN + if(u), 9{u + vr) = 2tvK + 6{u). (2.10) 



resulting from (2.2) and (2.4) 



Corollary 1. The set of distinct ordered in the sense (1.13) solutions of the trigonometric 
weak GBE (0.1) is in one-to-one correspondence with the equivalence classes T G Z M /7 



of integer vectors I G Z M factorised over the lattice of the periods 7 generated by the basis 
ln (n=l,...,M): 



M 



(ir. 



II 



r 



,n=l 



r g 7. 



(2.11) 



(2.12) 



Having thus characterized the distinct ordered solutions to the weak GBE in terms of 
the integer vectors J we proceed to characterize the solutions to the strong GBE. 

Definition 1. An equivalence class X G Z M /7 is called regular if all its representatives 
I lie off diagonals of Z A/ ; that is I m 7^ /„ for m ^ n. 

Theorem 3. The set of distinct ordered solutions to the strong GBE in the trigonometric 
case is in a one-to-one correspondence with the set of regular equivalence classes T. 
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Proof. The exclusion principle (proposition [2j demands to disregard the solutions v of 



the logBE (1.14) having multiple roots. More precisely, we discard the whole equivalence 



class V G !R." ra 7(7rZ) M if it contains a representative v G V such that v m = v n for some 
m 7^ n. Indeed, in such a case, by the monotonicity principle (proposition [3]), for the 
corresponding integer vector J we have I m = I n for the same m ^ n. Eliminating 
the equivalence classes generated by such representatives / we are left with the regular 
classes X. Vice versa, by the same monotonicity principle, if X is a regular class, all its 
representatives i" give rise to solutions v with distinct coordinates. ■ 

2.2 The rational case 

We now discuss the existence of solutions in the rational case. We will show that the 
condition of existence of solutions leads to a quite different description of the set of 
admissible integers J leading to distinct solutions of the rational GBE. 
We start by introducing the rational analog of the function x(w,c): 

u 

p(u,c)=iln( V ^)= ^^dv, c>0. (2.13) 
\u-icj J V 2 + cr 

— oo 

It satisfies p(u, c) + p(—u, c) = 2n and has asymptotics 

( s ( 0-2CU- 1 + 0{u- 3 ), u y oo , .* 

^' c ) = ( 2.-2^ + 0(^3), u^+oo ■ (2 ' 14) 

The integral of p is explicitly computable and reads 

-2c In \u\ + 0(1), u^—oo 

2-ku — 2clnu + O(l), u — Y +oo 

(2.15) 



\ pit, c) dt = ttu+2u arctan — cln I 1 H — - ) 
Jo c V c / 



In the rational case, <p and 6 are given by 

N K 

(p(u) = 2vr^ 00 + ^p( M + Re(en),Im(e„)) , 9(u) = ^ p(u + Re{ Vj ), Im(^)) . (2.16) 



n=l j=l 



Again 0(0) = nK follows from (1.2). We also remind that = e 2l7Tlfia 



Proposition 7. Assume that the phase (p^ is such that 

- 1 < (M + 1)^0 < 0. (2.17) 

then there exists a solution to the logarithmic GBE in the rational case if and only if for 
any subset J C {1, . . . , M} of cardinality \ J~\ = r with r = 1, . . . , M, the M integers 
(Ix, . . . , Im) satisfy the set of 2 M inequalities: 

T V~ l > K < J^ 1 ! ^ rMK- T ( T + l "> K + rN - 1. (2.18) 



11 



Corollary 2. (follows from the above result combined with propositions |l]|4j). The distinct 
ordered solutions to the rational generalized Bethe equations are in one-to-one correspon- 
dence with monotonously increasing sequences of M integers I\ < ■ ■ ■ < lu satisfying the 
set of 2 M constraints (2.18). 



Proof [Proposition [7]. Similarly to the proof of proposition |5j it is enough to prove 
the divergence of the potential W at infinity so as to ensure the existence of solutions. 
Reciprocally, if one is able to show that the potential is unbounded from below, then it 
cannot have a local minimum. Indeed due to the strict convexity of W, the local minimum 
would be global, contradicting the unboundedness of W. 

We first assume that the set of 2 M conditions (2.18) holds and prove that W goes to 



+oo along any ray. This is enough to prove that W admits at least one minimum. 

Next, we assume that at least one of the 2 M conditions (2.18) does not hold. Then we 



construct a ray along which the potential goes to — oo. As discussed above, this implies 
that W cannot have any local minimum. 

Suppose that the M integers Ii,...,I M satisfy to the constraints (2.18). It follow 



from the asymptotics (2.15) of the antiderivative of p(w, c), that $(u) and G have the 
asymptotics 



<S>(u) 



<p(t) dt 



0(«) 



9(t) dt 



2n<p 0O u - 2 (J2f =1 In | u | + 0(1), 

2u{ Voo + N)u-2 (j2f=i Im (0)) ln « + 0(1), 

-2 (Ef =1 Im(r/,))ln| M |+0(l), 
2nKu - 2 (J2? =1 M^)) lnM + 0(1), 



u — > — oo 
u — > +oo 



u — > — oo 
u — > +oo 



Up to a permutation of the 7,'s, W is a symmetric function of the v^s. Therefore we 



permute the coordinates of v with a £ & M in such a way that v a = (v^m, 



v, 



where t > and u : ^ i=1 u\ = 1 has its coordinates ordered in the following way 



<t{M)j 



tu, 



u 





with U\ < U2 < ■ ■ ■ < Ug < < U q+ \ < ■ • ■ < u s 



(2-19) 

Note that we have X p > for p = 1 . . . s and M — Ai + • • • + A s . The above permutation 
of the coordinates of v results in a permutation of the integers Ij h-> I^q). 
For further convenience we introduce the shorthand notations 



k 

p =i 



and Je 



j=Wt-l+l 



(2.20) 



and agree upon wo = 0. After some algebra, we recast W in the form 



W(v) = h®(tu k ) + - AiAj-e^ - t Uj ) + vr(M - l)K x ktu k - 2tt tu kJk 



k=i 



«J=1 



fe=l 



k=l 



(2.21) 
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Here, we can already send t — > +00 and compute the asymptotic behavior of all functions 
without problem: 



2^00^ + 0(111*), ke{l,...,q}, 
2^(^00 + N)tu k + 0(ln *) , ke{q + l,...,s}, 



(2.22) 



Xi\jQ(tUi — tuj) = 2nKt XiXj(v,j — Uj) + 0(ln*) . 

i,j=l l<i<ji'<s 

As a result, the asymptotics can be recast in the form 



(2.23) 



W(lu) = 2nt ( Y "i^i- ] + Y Vi T '-- 

k=q+l 



( + ) 



O(lnt) 



vfe=l 



where the two sequences 7^ read 



T k — Aft 



Jlc. 



K 

¥00 + y (A fc - 1) + ZfWfc-i 
^ + N + MiT - — (A fc + 1) - K(M - u;*) 



•A 



(2.24) 



The claim will follow as soon as we prove that the sum in the brackets in (2.24) is 



strictly positive. We first focus on the case where u q 7^ 0. Then the "(+)" and "(— )" 
partial sums are both strictly positive. Indeed we observe that Y^k=i u k T k~^ can be re-cast 
as: 



-(-) 



(2.25) 



k=i 



k=l 



0=1 



There, we have from the very definition of the ordering of the u k s that u k — u k+ i < and 



u q < 0. Also 



j=l p=l p=l j<p j=l 



K 



K 



u-'k 



Voo^k ~ —w k + — Y V\7 ~y^j< 



K , . K , 

< PooWk + —w k {w k - 1) - —w k (w k - 1) = y?ooU>fc < 



where we have used the Ihs of (2.18) and the definition of w k (2.20) so as to obtain the 
last line. It thus follows that Yll=i u k T k~" > > 0- 

It now remains to show that 2~2t=q+i u k T k + ^ > 0- Very similarly, we decompose 



s s 

^ / \ / \ !+) \ . 

k=q+l k=q+2 \j=k / j=<?+l 



Y UkT k +) = Y( Uk ~ Uk -^ J2 t j +) + u «+* ■ Yl 



T 



(+) 



(2.26) 



13 



As u k — Uk-i > and u q+ \ > 0, we ought to check the sign of the partial sums involving 
the Tfc . We set w k = Y^i=k = M — Wk-\ and get 



E 

j=k 



(+) 



K 

Wki^oo + N + MK)-—w k (w k + l)- Ia{j)> w k <foo + l>0 . (2.27) 

j=w k -i+l 



M 



The last inequality follows by applying the rhs part of (2.18). 

Now we assume that u q = 0. As m lies on the unit sphere it has to have non-zero 
coordinates. This means that necessarily s > 1. In its turn this ensures that even if the 
contribution of u q to the asymptotics is zero, there are still other ugs such that taken as 
a whole, the sum in the brackets in (2.24) is strictly positive. 

1, . . . , M satisfy the 2 M inequalities 



Therefore, provided that the M integers Ij, j 



(2.18), we see that the potential W(v) grows in all directions. This implies that it admits 
at least one local minimum, and hence that there exists a solution to the logarithmic 
Bethe equation (1.14) associated with these integers. 

Reciprocally, suppose that one of the 2 M inequalities (2.18) is not satisfied. We dis- 
tinguish between two cases. First suppose that there exists a subset J of cardinality 
r G {1, . . . , M} such that the Ihs part of inequality (2.18) is not satisfied. Then, as the 
inequalities only involve integers, we necessarily have 



K 



r(r — 1) 



(*) 



(2.28) 



fe=i 



Here a G &m is a permutation that maps {1, . 
vector v to infinity along the ray 

v a (k) = — 7= , for k = 1, . . . , r and v a ( k) 



, r} onto J . We choose to send the 
= , for k = r + l,...,M . (2.29) 



Then according to (2.24), we get 

W(v) 

where 



2nt 



.(-) 



O(lnt) 



.(-) 



K . 

r^oo + yr(r - 1) 



k=l 



> rif^ + 1 > . 



(2.30) 



(2.31) 



This shows that W is unbounded from below and hence cannot have any minimum. 
Similarly, if there exists a subset J of cardinality r G {1, . . . , M} such that the rhs of 
inequality (2.18) is not satisfied, then 

r 



r(r + 1) 
> rMK - V ' K + rN 



(2.32) 



k=l 



Here a G @m is any permutation that maps {1, 
vector v to infinity along the ray 

v a (k) = 4= , for fc = 1, . . . , r and v a{k) -■ 



. . , r} onto J . We choose to send the 
, for k = r + 1, . . . , M . (2.33) 
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Then according to (2.24), we get 

27Tt 



W{v) = ^=r 2 (+) + O(lnt) . (2.34) 



Paying attention to the conventions introduced in (2.19), we see that Ai = M — r, A2 = r 
and thus W2 = M, 

4 +) = r^oo + r(N + MK) - ^r{r + 1) - ^ I a(k) < r Voo < . (2.35) 

fc=i 

This shows that W goes to — 00 along that ray and the claim follows. ■ 

3 Number of independent solutions in the trigono- 
metric case 

3.1 Properties of factor modules 

Factor modules over Z M provide a convenient language for describing set of integers that 
are identified modulo shifts as those appearing in theorem |3j 

Let MatA/(Z) be the set of M x M matrices over the ring Z. Given a matrix G G 
Matjw(Z), Gm = "fu , Ike £ Z, one defines (G) to be the submodule in Z M generated by 
the rows of G. Then, the factor-module G = Z M / (G) is the set of equivalence classes 
of the integer vectors (Ji, . . . ,1m) differing by a vector from (G). In other words, G is 
generated by the generators forming a column vector e = (ei,...,e m )* subject to the 
relations 

M 

Ge = 0, or explicitly 7mne n = 0, m = 1, . . . , M . (3.1) 

n=l 

The factor module G relevant for counting the solutions to the generalized Bethe 
equations, is given by a matrix G with entries 

l a n = m 

Iran = S , ' / or "/ mn = a5 mn + 6(1 - 5 mn ) , (3.2) 

provided that one takes a = —K and b = N + (M — 1)K. 

Indeed, the above factor-module G is identical to the set of equivalence classes of the 



integer sequences / with respect to the transformations (2.9) 



The strategy we adopt for counting the number of solutions of the generalized Bethe 



equations (0.1) in the trigonometric case is as follows. We consider a factor module 



defined by the matrix G given in ( |3.2 ). We start by determining the total number 



of elements belonging to G. Then, in order to relate #G to the number of distinct 
solutions as described in theorem [3] we should discard all sequences I = Ylrk=i ^ e fe e ^ 
which have at least one repeating component Ij = I k , for some j 7^ k. Finally, as 



15 



the independent solutions of the generalized Bethe equations are in bijection with the 
monotonously growing sequences I\ < ■ ■ • < Im, we have to divide this last result by the 
number of allowed permutations so as to obtain the number of solutions. 

Recall that the permutation group acts in a natural way on the module Z M : 

a : (Ji, . . . , I M ) i-> (I<r(i), ■ ■ ■ , I<t{M))- (3.3) 



It is obvious that the submodule (G) given by the matrix G (3.2) is invariant under the 
above action of &m- The induced action of &m can be then defined on the factor-module 
G. 

An orbit of a point I = J2jLi Ij e j G G is the set of all I CT = J2fLi ^{j) e j f° r ° e ©m- 
A regular, or generic orbit is an orbit that has the maximal possible number of distinct 
points, that is Ml. An example of points I = Ylk=i h e j no ^ giving r i se to regular orbits 
are the so-called diagonal points, that is points having at least two integers coordinates 
equal: Ij = Ik for j ^ k. The regular orbits of the action of &m on G are described by 
the following lemma 

Lemma 1. The regular orbits of the factor module G subordinate to the period matrix 
G (3.2) are generated by all elements I = X^fcii Ik e k such that the integers Ij form 
monotonous sequences I\ < ■■■ < Im and the vector (Ii, . . . , Im) belongs to the fun- 
damental domain & of G 

& = l^x e R M : Xi = U(a - b) + b^tp, ^e(0,l)|. (3.4) 

Proof. It is clear from the previous discussion that the regular orbits can only be given by 
off-diagonal elements. Moreover, one can always represent G by points in its fundamental 
domain. Let I = Y^h=i be defined in terms of a monotonous sequence of integers 
Ji < I2 < • • • < Im such that (I±, . . . , Im) lies in the fundamental domain of G. 

Assume that the ©M-orbit of I in G = Z / (G) is not regular. Then, there exist at 
least two distinct permutations a and a' such that I a = 1^, where I CT = Y2 P =i 

I a{p) e p . In 

other words, 



sM 

'1 



M 

(J CT (1), . . . , Icr(M) ) = {Ia'(l), • • • , Ia'{M)) + %(7jl> • • • > 7jAf) , (3.5) 



for some a,- e Z and 7 mn given by (3.2). Applying the inverse permutation, we get 



M 



(It, . . . ,I M ) = (J CT ' CT -i(l), • • • , Ia'a-i-iM)) + 2j a <T- 1 (j)('y<T- 1 (j)a- 1 (l), ■ ■ ■ ,1^(J)(T- 1 (M)) ■ 

3=1 

(3.6) 

The period matrix G (3.2) can be recast as 

G = (a - 6)id + ■ y*, 
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where x t — y* = (1, . . . , 1) are permutation-invariant vectors. Clearly, this implies that 
the period matrix is invariant under permutations, ie j a (j)o-(k) — Ijk- It thus follows that 
there exist integers (3j G Z and a non-trivial permutation n G &m such that 



M 



I P = h( P ) + (a - b)P P + bJ2&- (3-7) 



Summing (3.7) over all the p's we get that 

M M 

(a + (M-l)b)J2fa = ° X^ = ' ( 3 - 8 ) 

£=1 £=1 

the latter stemming from the fact that the period matrix G is assumed to have a non- 
vanishing determinant, cf (3.17). Thus, we get an expression for j3 p in terms of the 
parameters t p describing the points I p , namely 

0p — t p — t 7r ( p ) G (—1,1). (3.9) 

As P G Z, we necessarily have f3 p = for all p. As I is an off-diagonal element 
(Ii,...,Im) = (^Tr(i)) • • ■ j In{M)) implies that n = id, which contradicts the fact that 
a ^ a'. m 

The above lemma leads to a complete description of the solutions of the generalized 
Bethe equations in the trigonometric case. 



Proposition 8. The solutions of the trigonometric generalized Bethe equations (0.1) are 
in a one-to-one correspondence with the regular orbits of the factor-module 

G = TL M j (G) 

with respect to the action of permutation group &m- 

In the remaining part of this Section, we will prove 



Theorem 4. For a matrix G of the form (3.2), where a and b are not necessarily given 
by (2.7), the number of regular orbits of the factor module G = 7L M / (G) is 

, [0)6] \a+(M-l)b\ (\a-b\-l 



m " K m-i ) (3 ' 10) 

Taking a = N + (M — 1)K and b = -K we get that 
Corollary 3. The number of solutions to the generalized Bethe ansatz equation in the 



trigonometric case is given, for N = 1, by the Fuss-Catalan numbers , and by 

i(k,n) N [KM + N - l\ N [KM + N 



° M M\ M-l ) KM + N\ M J' ( 3 ' n ) 

for arbitrary N. Note that if C(x) = YI'mZo^m' 1 ^^ 1 ^ s the generating function of the 
Fuss-Catalan numbers then C N (x) = Y^mIo^m^'^^ is the generating function of the 



numbers C M 



(K,N) 
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3.2 Number of elements in G 



Proposition 9. Let G be a factor module of Z M defined by some period matrix G G 
Matjvf(Z) such that detG ^ 0, then G is a finite set, and its cardinality, up to an eventual 
sign, equals to the determinant of G: 



#G = | detG | 



(3.12) 



Proof. It follows from elementary algebra that the matrix G admits a Smith normal 
form: 

G = UDV, U,D,V G Mat M (Z) , D = diag(di, . . . , d M ) . (3.13) 
and the matrices U and V are invertible (it det U = ±1 , det V = ±1). 



Substituting (3.13) into the set of relations (3.1) determining G we obtain 

UDVe = . 



(3.14) 



The matrix U being invertible, we multiply the relations by U^ 1 from the left and 
obtain an equivalent system of relations. On the other hand, due to the invertibility of V, 
the system of generators e for the module Z M can be replaced by the equivalent system 
of generators e' = Ve. 

As a result, we get an equivalent description of the factor- module G in terms of the 
generators e' and the relations 



De! = 



ie d m e m 



0, for m 



M 



(3.15) 



Thus, the factor-module G decomposes into a direct sum of rank 1 modules, each with 
the single generator e' m and the single relation d m e' m = 0. Obviously, each component 
consists of \d m \ elements that is equivalence classes whose representatives can be chosen, 
for example, as 



0, Cm; 26 m ; 3Cm; ■ ■ ■ j (|<^m| l)*3r. 



Thus the cardinality of G equals Y[ m \d r , 
IdetDl. 



(3.16) 

[det D\. It remains to note that, by virtue 



of (|3_13j) |detG| = 
Corollary 4. Let 



be the factor module defined by the period matrix G (3.2). Then 
#G= |detG| = | (a + (M — 1)6) (a — 6) M_1 1 . 



(3.17) 



Proof. In the case of the matrix G of interest ( |3.2 ), detG can be computed thanks to 
the well-known formula for the determinant of rank-1 perturbation of a matrix. Let A be 
an arbitrary square matrix, x be a vector-column, and y l be a vector-row. Then 



det (A + xy l ) = det A + y*Ax, 

where A is the adjoint matrix that is AA = det A. 
In the case of interest we have G = A + xy*, where 

A = (a - 6) idjur, x= j , y* = 

w 



(3.18) 



(3.19) 
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As we have already pointed out, formula (3.17), constitutes only the first step towards 



the answer; what we need to count are not all the points of G but the regular orbits in 
respect to the action of the group &m- 

3.3 The diagonal submodules of G 

In order to describe the diagonal submodules of G, we start by presenting an example. 
Example. When M = 2, the period matrix 

G =(lS < 3 - 20 > 



produces the system of relations 



ae 1 + be 2 = (3.21) 
oei + ae 2 = (J . 



Diagonals correspond to I\ — I 2 = ^12 that is to the points /^(ei + e 2 ). In other words, 
the element e\ + e 2 generates a submodule, which we shall call a diagonal submodule. The 



relation for this submodule is obtained by adding together the equations (3.21): 

(a + 6)(ei + e 2 ) = 0. (3.22) 

The size of the submodule (the number of points on the diagonal) is obviously \a + b\. 

Let V be a partition of the set {1,. . . , M} into J disjoint non-empty subsets V = 
{Vi, . . . ,Vj}. To each partition V we associate the diagonal submodule V-p generated 
by J generators e-p. = Yliev- e «- ^ s ma trix of relations is obtained from the matrix G by 
adding the rows labeled by i G Vj. 

Example. For M = 3 there are 5 partitions of the set {1,2,3}: the 3-component one: 
{{1}, {2}, {3}} corresponding to the original module (G), and 4 other partitions generat- 
ing diagonal modules: three 2-component ones: {{1,2}, {3}}, {{1,3}, {2}}, {{2,3}, {1}} 
and one 1-component: {{1,2,3}}. 

The 2-component partition {{1, 2}, {3}} gives rise to the generators {ei + e 2 = ei 2 , 63} 
and relations (a+6)ei 2 + 26e3 = 6ei 2 +ae3 = 0, similarly for the rest two. The 1-component 
partition {{1, 2, 3}} gives rise to a single generator e\ + e 2 + = ei 2 3 and a single relation 
(a + 26)ei 23 • 



Proposition 10. LetV = {Vi, . . . ,Vj} be a partition of {1, . . . , M}. Then the associated 
diagonal submodule D-p is isomorphic to the diagonal submodule T>\ corresponding to the 
partition V\ = {1, . . . , Ai}, . . . , Vj — {M — Aj + 1, . . . , M}, where Xj = #Vj gives the 
number of elements in the partition Vj . 

Proof. Taking the cardinalities j^Vj of the components of the partition V of the set 
{1,...,M} and rearranging them in the decreasing order we obtain a partition A = 
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{Ai > A2 > ■ ■ ■ > Aj} of the number M into J nonzero numbers: | A| = Ai + . . . + Aj = M. 
Clearly, partitions V corresponding to the same partition A are related by permutations 
of the numbers {!,..., M}, and the corresponding diagonal submodules are isomorphic. 



Corollary 5. The number of elements in every diagonal submodule T>-p only depends on 
the number J of components of the partition V = {V±, . . . , Vj} and is given by 



#V P = Af — (a + (M — 1)6) (a - b) 



j-i 



(3.23) 



Proof. As T>v and T>\ are isomorphic, is is enough to count the number of elements if 
T>\. The matrix elements of period matrix G\ for the diagonal submodule T>\ are given 
by 

9% = W(° + (A, - 1)6) + (1 - 6 jk ) A,6, j, k = 1, . . . , J. (3.24) 



Rewriting (3.24) as 

g% = {a-b) S jk + \jb, (3.25) 

we notice that G\ decomposes into a sum of a diagonal matrix A = (a — 6) idj and a rank 
1 matrix xy*, where Xj = A,- and y k = 6. Applying equation (3.18) we obtain the claim. ■ 



3.4 Elimination of diagonals and Stirling numbers: proof of the- 
orem m 

We are now in position to prove Theorem |4| 

The case M = 2 was analyzed at the beginning of the previous Section. There, one 
has the decomposition 

2!Cf' b] = |detGn| - |detG 2 | = |A^| - |A?| = \a + b\ (\a-b\ - 1). (3.26) 

For M = 3 we have 3 partitions of the set {1, 2, 3} of the type A = (21). Subtracting 
the corresponding points count from the number of points in the original module G we 
obtain 

|det Cm I — 3 |det G%i\ ■ 

Note, however, that each of the 3 diagonals I\ = I2, I\ = I3, I2 = h contains the 
common sub-diagonal I\ = 1% = I3. We have thus subtracted it thrice. To compensate, 
we need to add it back twice. Finally, we get: 

3!Cf' b] = IdetGml -3|detG 2 i| + 2|detG 3 | 
= |Ajj| - 3 |A|| + 2 |A?| 

= \a + 2b\ (\a -b\- l)(|a - b\ - 2). (3.27) 
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Proceeding similarly for M > 3 and using corollary [5] stating that det G-p only depends 
on J we arrive to the following expansion: 



M 



M\C^ b] = J2(-V M - J SMj\A 



M I 
J 



J=l 



M 



+ (M -l) b \J2(-l) M - J SMj\(a-b) 



j-ii 



(3.28) 



j=\ 



where smj are some combinatorial coefficients that do not depend on a and b. 

To determine smj we invoke their independence on a and b and set b = 0. Then, the 
factor-module G becomes a hypercube having sides of length \a\. The regular orbits are 
labeled by the sequences < Ji < 1% < ■ ■ ■ < Im < M — 1- The number of such points is 
obviously 

.[0,0] = f\a\\ = H(|a| - l)...(|a| -M+l) 
ImJ M! 



c 



(3.29) 



Setting b = in (3.28) and substituting (3.29), we obtain its expansion into powers of 



(H-i)...(H-m + i) = 5>i; 



M-J„ \ n \J 

smj \a\ 



(3.30) 



j=i 



The latter coincides with the definition of the (unsigned) Stirling numbers of 1st kind 
PHQ2J. In Knuth's notation pi]: 



M 
J 



(3.31) 



Returning to the general case (3.28) and using the definition of the Stirling numbers 
we obtain 



M 



M\C 



[a,b] 
M 



+ (m-i)&i5>i; 



M-J 



J=l 



M 
J 



I (0-6) 



j-11 



\a + (M — l)b\ 



\a — b\(\a — b\ — 1) . . . (|o — 6| — M+l) 



|(a-6)| 

\a + (M - l)b\ {\a - b\ - 1) . . . (\a — b\—M+ 1). 



The last expression coincides with (3.10). 



4 Proof of theorem 2 



In this Section we prove theorem [2} The proof is carried out by establishing a bijection 
between the solutions to the generalized Bethe equation in the trigonometric and rational 
cases. 

One might expect that the solutions to the generalized Bethe equations in the trigono- 
metric case with parameters {er] k } and {e^}, scale as 0(e) in the e — > + limit (up to 
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some possible shifts by (7rZ) M ). If this is indeed the case, then the rational equations 
can be considered as a degeneration of the trigonometric ones. Of course, one should be 
rather careful in such a limiting procedure. When e — > + , it could be that some solu- 
tions to the trigonometric equations may not go to zero, and thus one might loose some 
solutions in the procedure. In particular, there might not be a one-to-one correspondence 
between the two sets of solutions. We prove that this does not happen: one can build, 
for e small enough, a one-to-one correspondence between the solutions to the generalized 
Bethe equations in the trigonometric and rational cases. 

From now on, we assume that the parameters {^} and {f]k} occurring in the trigono- 
metric case are rescaled by e > 0. Therefore, the trigonometric functions F and S are 
given by 

F(u) = F 00 T[ Sin / (M + £ l n \ where F^ = e ~ 2inLPc ° , -1< (M + l) Voo < 

s(u) = f\ sin{ ; +£r]n \ . 

In the following, the parameters {rjk}^ =1 and {^ r }^ =1 are to be considered fixed. The 
only varying parameter will be e. We will call it the dilatation parameter. From now 
on, we shall only focus on the solutions v to the trigonometric GBE belonging to the 
fundamental domain (— n/2, tc/2], ie for all j € {1,...,M}, Vj G (— 7r/2, 7r/2]. We also 
introduce the following definitions for the function (p £ and 9 £ : 

N 

ip £ {u) = 271-^00 + ^2x(u + £ Re(Cn), £ Im(O) > 

n=l 

K 

°e(u) =^2x(u + eRe(r] j ),£lm.(r} :j )) , 



where x{ u -> c ) ls given by (2.1). 



4.1 Rational constraints on the trigonometric integers 

In this section we prove the following theorem 

Theorem 5. There exists Eq small enough such that any sequence Ii, . . . , Im defined by a 



solution to the trigonometric GBAE (1.14) with dilatation parameter e e (0,e ) satisfies 
the inequalities: 

Nr-l- K r{T ^ ^ + KMr >^h> # ~ , Vjc {1.....M}. (4.1) 

The main consequence of this theorem is that any solutions of the trigonometric GBE 
with sufficiently small inhomogeneities {erjk}, gives rise to a set of M integers (Ii, . . . , Im) 
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satisfying the set of inequalities issued from the minimum condition of the rational po- 
tential. This means that there is at most as much solutions to the trigonometric GBE as 
there are solutions to the rational GBE. 

To prove the theorem we need to establish two preparatory lemmas. We first observe 
that for e small enough, any coordinates Vj of a solution v to the trigonometric logarithmic 
Bethe equations stays uniformly away from the boundaries ±7r/2. 

Lemma 2. There exists an Eq > and an a > such that given any solution v = 
(vi, . . . ,%) to the trigonometric logarithmic GBE (1.14) with dilatation parameter e G 
(0, e ), one has 

for j 1 M . (4.2) 



Vj e 



IT TT 

— h a, a 

2 ' 2 



Proof. We prove the statement by contradiction. Assume that there exists a sequence 
e n — > and a sequence of solutions to the logarithmic GBE arising in the trigono- 
metric case and corresponding to a dilatation parameter e n such that 

or v\,> — > - (4.3) 

n— >+oo z n— >+oo Z 

Here, we choose to work with solutions having ordered coordinates —tt/2 < < 
■■■ < v$ < tt/2. As the sequence is bounded in norm, it admits a converging 
subsequence. We continue to denote this converging subsequence by v^ n \ It fulfills 

vi n) — > -tt/2, se {!,...£} 

vi n) — > G (-tt/2, tt/2), se{£ + l,...M-£'} (4.4) 

— ► tt/2, s G {M - £' + 1, . . . , M} . 

n— >+oo 

Due to (4.3), at least one of the Vj has to converge to one of the boundaries points ±7r/2. 

Hence, £ + £' > 1. Taking the product of the trigonometric GBE for where s runs 
through the set S = {1, £} U {M - £' + 1, . . . , M} yields 

M M 

i = n f ^ ■ n n - ^ = n ^ ■ n n - ^ n) ) • ( 4 - 5 ) 

ses seSfc=i ses sG5fe=i 

In the last equality we made use of the space parity condition satisfied by S so as to 
cancel out the products S(vs — v k )^( v h ~ f° r k,s £ S. It is easy to see that for 
s G S, F(vi nS) ) — > Fqo. As for k S, converges to some point in (— tt/2, tt/2). We 
have, 

v^-v^ — > 8v sk G (-tt,0) U (0,tt) Vs G <S and jfe G {1, . . . , M} \ <S . (4.6) 



n— >+oo 



In its turn, this means that S(vi n ^ — v^ 1 ) — > 1. Thus, the n — > +oo limit of (4.5) 
yields 
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As 1 < £ + £' < M + 1, this contradicts the fact that -1 < (M + l)y?oo < 0. ■ 

We now observe that, if we take any two points v,w lying in (— tt/2 + a, tt/2 — a) 
for some a > 0, then 9 e will be almost contained in the interval [0, 27r], ie 8 e (v — w) G 
[—u, 2ir + v], where v — )■ as e — > 0. 



Lemma 3. Let < a < tt/A, then there exists Eq > such that for all e G (0,£q) 



2irK - 2tt^ > 9Jv - w) > 2tt^°' 



M 2 



M 2 



V v , u> G 



7T 7T 

ha. a 

2 ' 2 



(4.8) 



Here, we remind that — 1 < < 0. 
Proof. Let Imf? > 0. Then, as \ ls strictly increasing on the real axis, given v,w G 
[— 7r/2 + a, tt/2 — a] we have 

x(— n + 2a + e Re T], elm rj) < x( v ~~ w + £Rerj,elmr]) < x( n ~ 2a + eRer],elmrj) . 

In the e — > + limit, the //is goes to zero whereas the rhs goes to 2tt. Indeed, taking e 
small enough so that < a — e\ Re(r/)|, one gets 



x(— vr + 2a + e Re(?y), e Im(r/)) 



-7r+2a+eRe(??) 



sinh(2e Im(?7)) 



< 



sin(w + z£ Im(?y)) sin(u — z£ Im(^)) 
sinh(2e Im(r/)) ■ 7r/2 



sin(a + 2£ Im(?7)) sin(a — 1111(77)) e ^ 0+ 
Similarly, using the quasi-periodicity of x we get 



du 



. 



x(tt — 2a + eRe(?7),5lm(r7)) — 27r 

-2«+eRe(rj) 



sinh(2e Im(?])) 



du 



sin(w + ie Im(r/)) sin(w — ie Im(r/)) e^o+ 



It thus follows that there exists an e > such that for all 

k G {1, . . . , i^} and v,w G 

we have 



7T 7T 

— ha, a 

2 ' 2 



< x(f - u> + £Re(7/fe),£lm(7/ fc )J < 27r 



(4.9) 
(4.10) 



KM 2 ^ v v ' n v ' y/ KM 2 ' 

uniformly in e G (0, £0). Summing these inequalities over k we get the claim. ■ 
We are now in position to prove the main result of this subsection: 

Proof [ Theorem [5] . We choose an £0 > such that the conclusion of lemmas [2] and [3] hold 
simultaneously. Let e be such that £0 > £ > and v be any solution of the trigonometric 
logarithmic GBE with parameters {er] k } and {e£ p }. Then, according to the conclusions 
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of lemma [2] there exists an a such that Vj E (—ix/2 + a, n/2 — a). By virtue of lemma [3j 
we have 

27r v^ < d e(vj - v k ) < 2nK - 2vr^ . (4.11) 
M 2 M 2 

Also, since Vj G (— 7r/2 + a, 7r/2 — a) and one can always assume e to be small enough 
so that a — e\ Re(^)| > for k = 1, . . . , N, we have 

27r(iV + ^oo) > ^ e (^) > 271-^00 . (4.12) 

Let J" C {1, ... , M} such that # J" = r, then 

2tt h = + Yl Yl *( v * ~ v i) 

kej kej kej j^k 

M 

r(r — 1) r(M — r) 

> 2ttK V 2 ; + 2ixr Voo + 2tt V M2 ; ^ 

> 27r(r + l) Voo + 27rAr r ^~ ^ . 

It remains to use the bound — 1 < (r + l)y?oo < and invoke the fact that ^2 k&s h is 



an integer so as to obtain the equality in the rhs of (4.1) 



00 



One gets the second set of inequalities in a similar way 

27T 4 = Y [ 6 e( V j - V k) + 9e(Vk ~ Vj)] + Y V^k) + 6 ^ Vk ~ V ^ 

kej j,keJ kej kej j^J 

j>k 

< 2nK— + r2ir(p 00 + 2nNr + 2irK(M - r)r - 2tt^— — — tp 

2 M 2 

< 2n (^Nr - 1 - y ^ + KMr^j . 

When passing to the last line, we again have used the fact that Ylikej ^ k ls an integer. 



4.2 From the rational to the trigonometric solutions 

One can visualize the generalized Bethe equation as the set of constraints defining the 
real valued zeroes v with pairwise distinct coordinates of the mapping 

y {vm} :C M ^C M 

Z {Q{vUO( Z l I Z )> • • • > Q{riUt}( Z M I z)) . 

We remind that the trigonometric Baxter polynomial Q,{q\te\ ( z m I z ) has been intro- 



duced in (1.7). 
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Proposition 11. Let v be a solution of the rational generalized Bethe equations with 
parameters {r/k} and Then, there exists an e v > and a mapping g v (e) from the 

interval (—e v ,e v ) to some open neighborhood of in K M such that eg v (e) is a solution of 
the trigonometric generalized Bethe equations with parameters {en} and {e£}. 



Proof. The function Q given by 

£:(-!;!) x 

is continuously differentiable and 



■>M 



M 



(e, z) i — y is 



-N-MK 



N 



M K 



lim — • Q k (e, z) = <{ JjK + £ r ) • JJ JJ(^fc - v m + r] e ) 



d 2 



r=l 



m=l l=\ 



dzjdzk 



W(z) 



(4.13) 

Here W is the strictly convex potential (1.18) appearing in the rational case. As any 
solution to the generalized Bethe equations is real, the product in the pre-factor in the 
rhs is non-zero. We have shown in proposition [4] that W has a positive defined Hessian 

d 2 



det 



dzjdzk 



W(z) 



>0, 



(4.14) 



therefore detjy/ [d Zk Ge(0,v)} ^ and we can apply the implicit function theorem to the 
mapping Q in the vicinity of the point (0, v). 

There exists open neighborhoods (— £ v ,£ v ) of 0, and V v of v and a C 1 mapping g v : 
(—£„,£„) — )■ Ki, such that g v (0) = v and 

{(e,z) G (-£„,£„) x K : ^ {er?}){£S} (£2;) = 0} = { (e,g v (e)) : £ G (-£„,£„) } . (4.15) 

g v (e) is thus a solution of the trigonometric GBE with parameters {erj} and {££}• ■ 



Proposition 12. Lei t> and v' be two distinct solutions of the rational GBE. Let £q = 
min(£„, e v i) > 0, where e v is as given in the above proposition. Then, for any £ G (— £o, £o), 
9v{e) ^ 9v{e). 

Proof. As e i— >■ g v {£) is continuous, it follows that 

¥e {e\9v]k(e)) an d £ {e\g v ] k (e) - e\g v ] h (e)) (4.16) 
are all continuous in £ on (— £ ,£o). Hence, the combination 

M 

£^ Ve (e\g v ] k (e)) +^2e e (e\g v ] k (e) - e\g v ] k (e)) (4.17) 

is also continuous. However, for any £ small enough, eg v (e) solves the trigonometric 
generalized Bethe equations. Hence, the function defined in (4.17) is integer valued. Being 



26 



a continuous function of e, we deduce that it is constant. The value of this constant can, 
for instance, be determined by setting e — 0. It thus follows that given two distinct 
solutions v and v' of the rational GBE, one obtains two solutions of the trigonometric 
GBE that are characterized by two distinct set of integers I and I'. Therefore these 
solutions cannot coincide. ■ 



Hence, for < e < Eo the family {e g v (e)}, where v solves the rational GBE, provides 
as many distinct solutions of the trigonometric GBE with parameters {er}} and as 
there are solutions to the rational GBE with parameters {rj} and {£}. We have just 
established that there are at least as many solutions to the trigonometric GBE as there 
are to the rational ones. This ends the proof of theorem [2j 



5 Conclusion 

In this article, we have studied a generalization of the Bethe equation where the ^-matrix 
has several poles and zeroes in the complex plane. We have characterized the set of 
solutions in the so-called repulsive regime. In particular, we have provided a thorough 
count of the number of solutions to these equations. We showed that this number is 
given or closely related to the Fuss-Catalan numbers. This allowed us to establish two 
yet-unknown combinatorial interpretations of these numbers. On the one hand, we have 
shown that the Fuss-Catalan numbers count the number of integers satisfying certain 
inequalities and on the other hand they give the number of regular orbits in certain factor 
modules. 

For the moment, such generalized Bethe equation do not correspond to any integrable 
model. We plan to investigate integrable models giving rise to such generalized Bethe 
equations in the future. It also seems extremely appealing to analyse the combinatorial 



identities arising from other classes of Bethe equations then those given in (0.1). For 
instance, estimating the number of solutions to the multi-pole generalizations of Gaudin- 
type equations should also lead to combinatorial problems for counting integers subject 
to sets of constraints. 
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